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ABSTRACT 



Aims. M87 is a key object whose study can reveal the complex phenomena in cooling cores. We use a deep XMM-Newton observation 
of M87 to produce detailed temperature, pressure and entropy maps in order to analyze the physical processes of cooling cores and of 
their heating mechanisms. 

Methods. We employed both broad-band fitting and full spectroscopical one-temperature model analysis to derive temperature and 
surface brightness maps, from which the pseudo-deprojected entropy and pressure were calculated. We discuss possible physical 
interpretations of small deviations from radial and elliptical symmetry in these maps. 

Results. The most prominent features observed are the E and SW X-ray arms that coincide with powerful radio lobes, a weak shock 
at a radius of 3', an overall ellipticity in the pressure map and a NW/SE asymmetry in the entropy map which we associate with the 
motion of the galaxy towards the NW. For the first time we find evidence that cold, metal-rich gas is being transported out of the 
center, possibly through bubble-induced mixing. Several edges in the abundance map indicate an oscillation of the galaxy along the 
NW/SE direction. Furthermore, the radio lobes appear to rise along the short axis of the elliptical pressure distribution, following the 
steepest gradient of the gravitational potential, and seem to contain a nonthermal pressure component. 

Key words. galaxies:individual:M87 - galaxies: intergalactic medium - cooling flows - X-rays:galaxies:clusters 



1. Introduction 

As early as 1973, the first detailed X-ray observations with the 
Uhuru and Copernicus satellites revealed that the X-ray luminos- 
ity at the center of some galaxy clusters and elliptical galaxies 
has a very sharp peak, implying a high enough centra l gas den- 
sity th at the cooling time falls below the Hubble time dLea et al.l 
1973). Consequently, the so-called cooling flow model was de- 
veloped which asserted that the gas, in losing energy by means 
of radiation and in the absence of other heat sources, cools down 
to temperatures so low as to become undetectable in X-rays , 
and sinks towards the cl uster core (e.g. lFabian & Nulsenl (1 19771) . 
ICowi e & BinneyJ (119771) ). The model however failed to explain 
the ultimate fate of the cooling gas, and the evidence in wave- 
bands other than X-rays showed that much less gas was condens- 
ing i nto a cooler state than the "coo ling flow" model predicted 
(e.g. lMcNamara & Q'Connelil(ll989l) ). 

The latest generation of X-ray satellites, in particular XMM- 
Newton and Chandra, provided more detailed spectral informa- 
tion showing that indeed much less gas undergoes such a "cool- 
ing flow" than the initial model based on the X-ray luminos - 
ity peak predicted dPeterson et all (l200lh |Peterson et al.l (120031) . 
IBohringer etaf] (l200ll) .lMatsushit a etal.1 d2002l) ). This led to the 
conclusion that a heating source must exist which interacts with 
the gas in a fine-tuned way, stopping the coo ling at just the 
right l evel to reproduce the observations (e.g. Bohringe r et al.l 
(120021) '). Active galactic nuclei (AGN) are found in most clus- 
ters previously believed to host a "cooling flow" (lBurnsl fl990) 
and represent the most favoured candidates for providing such 
heat input. The main challenge remains the understanding of 
the exact mechanisms by which the heating takes place (e.g. 



David et al.l d200ll).IChurazov et alJ j200lh . lFabian et ail d2003l). 
Brug gen et all d2005l) . lHeinz et al.l i20m ). 
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M87, the central dominant galaxy in the nearby Virgo clus- 
ter, is an ideal candidate whose study may shed light on the com- 
plex phenomena in cooling cores. M87 was recognized early on 
as having a peaked X-ray surface brightness corresponding, in 
the absence of heating, t o a cooling flow wi th a mass deposi- 
tion rate of 10 M /year dStewart et alJ 1 19840 . However, using 
the sp ectral information from XMM-Newton, IBohringer et al.l 

(2001) put an upper limit on the mass depositio n rate at least one 
order o f magnitude bel ow this value (see also [Matsushita et al.l 

(2002) , Molendi (2002)). The energy input which may prevent 
the cooling of the gas could be supplied by the central AGN, 
powered by the galaxy's cent ral supermassive bl ack hole which 
has a mass of 3.2 x 10 9 M Q dHarms et alj|1994l) . Based on the 
observed complex system of radio lobes, be lieved to be associ- 
ated with the AGN jet and unseen counteriet. lOwen et alJ {2000) 
showed that the mechanical power input from this AGN is more 
than sufficient to compensate for the energy loss through X-ray 
radiation. 

The radio plasma in M87 clearly interacts with the X-ray gas. 
The best proof of this is the fact that the two most prominent 
radio lobes extending to the east and south-west of the M87 cen- 
ter correspond to regions of enhanced X-ray surface brightness, 
also known as the E and SW X-ray arms. Moreover, these E and 
SW arms show a mo re complex multi-phase temperature struc- 
ture (Molendi 2002) while a previous observation with XMM- 
Newton showed that the rest of the ICM is loca lly well approx- 
imated as singl e -phase (Matsushita et al. 2002). This motivated 
IChurazov et al.l d2001l) to develop a model according to which 
the bubbles seen in the radio map buoyantly rise through the hot 
plasma uplifting cooler gas from the central region. 
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Fig. 1. M87 surface brightness map, from all three EPIC detec- 
tors vignetting corrected and combined, in the 0.5-7.5 keV en- 
ergy range 



Further fe atures in the X-ray su rface brightness map of M87 
described by For man et alj (|2005) from Chandra observations 
include a ring of enhanced emission with a radius of 14 kpc, 
possibly associated with a weak shock which may be one of 
the mechanisms contributing to the heating of the ICM, addi- 
tional arc-like features at 17 and 37 kpc, and excess emission 
regions NW and SE of the core. Using deeper Chandra data, 
Forman et al. (2006) confirmed the presence of a weak shock at 
14 kpc and determined its Mach number as M~1.2. 

Detailed temperature, entropy and pressure maps are a pre- 
requisite for the physical understanding of the phenomena be- 
hind these features. This is the motivation for the second, deeper 
observation of M87 with XMM-Newton, the results of which are 
presented in this paper. 

The paper is laid out as follows. In Section 2 we present the 
data set and data analysis methods. The features seen in the re- 
sulting temperature, entropy and pressure maps are described in 
Section 3. In Section 4 we present more quantitative analyses of 
these features and possible physical explanations. Our conclu- 
sions are summarized in Section 5. We adopt a reds hift for M87 
of z= 0. 00436 and a luminosity distance of 16 Mpc (iTonrv et al.l 
2001), which yields a scale of 4.65 kpc per arcminute. 

2. Observations and data reduction 

M87 was first observed with XMM-Newton in the performance 
verification (PV) phase on June 19, 2000, for 60 kiloseconds 
(ksec). A subsequent 109 ksec observation was performed on 
January 10th, 2005. We will focus in this work primarily on the 
EPIC data from the second observation, in which the PN detector 
was operated in the extended full frame window mode while the 
full frame mode was employed for the MOS detector. The results 
from the RGS data analysis for the two combined observations 
can be found in I Werner et al.l d2006l) . 

We extracted a lightcurve for each of the three detectors sep- 
arately and excluded the time periods in the observation when 
the count rate deviated from the m ean by more than 3<x in order 
to remove flaring from soft protons (iPratt & Arnaudl2002l) . After 



this cleaning, the net effective exposure is ~62 ks, ~79.3 ks and 
~ 80.5 ks for pn, MOS1 and MOS2, recording in total 13.2, 4.7 
and 4.9 million photons respectively. For data reduction we used 
the 6.5.0 version of the XMM-Newton Science Analysis System 
(SAS); th e standard ana l ysis m ethods using this software are de- 
scribed in Wat son et al.l d200ll) . A vignetting-corrected flux map 
of the observation for all three detectors combined in the energy 
band 0.5-7.5 keV is shown in FigureQ] 

For the background subtraction, owing to the extended na- 
ture of our source, which covers the entire field of view, we 
used a collection o f blank-sky maps from which point sources 
have been excised (Re ad & Ponmanl 120 03). We transposed the 
blank sky-maps to a position in the sky corresponding to the ori- 
entation of XMM-Newton during our observation of M87. The 
backgrounds were scaled according to the corresponding expo- 
sure times for each detector, after having removed the time peri- 
ods affected by proton flares. The scaling factors in this case are 
calculated as the ratios between the observation exposure times 
and the background exposure times. An alternative method for 
calculating background scaling factors was employed by requir- 
ing that in the higher energy region (10-12 keV for MOS, 12-14 
keV for PN) the emission in the outskirts of the M87 observation 
should have the same level as the background, since we expect 
little or no such high-energy emission from the gas halo around 
M87. The background scaling factors were in this second case 
calculated as the ratios between the total number of observa- 
tion to background counts in the higher energy region specified 
above, in a 7'-9' ring. The values for the scaling factors calcu- 
lated from the two methods agree to within less than 5%. This 
justifies confidence in having chosen a correct scaling. 

Where a stable minimum signal-to-noise ratio was needed 
for our analysis, in particular for determining the temperature, 
we employed an adapti ve binning method based on weighted 
Voronoi tessellations (iDiehl & Sta tler 2006) , which is a general - 
ization of the algorithm presented in lCappellari & Copinl d2003l) . 
The advantage of this algorithm is that it produces smoothly 
varying binning shapes that are geometrically unbiased and do 
not introduce artificially-looking structures, as can be the case 
for e.g. rectangular bins. Binning the dataset from the second 
XMM-Newton observation alone to 10 4 counts per spatial bin 
in the 0.5-7.5 keV energy range, which should yield tempera- 
ture values accurate to less than 5% for one-temperature model 
spectral fits, we obtain bins roughly corresponding in size to the 
extent of the XMM point-spread function (PSF). Therefore, in 
this paper we only discuss the results from the second observa- 
tional dataset, which contains enough counts for the purpose of 
the current analysis. More detailed spectral models and spectral 
fitting which are the scope of future work will use data from both 
observations. 

Out-of-time events were subtracted from the PN data using 
the standard SAS prescription. 

2. 1 . Broad-band fitting 

In order to obtain a first impression of the temperature and abun- 
dance distribution of the hot gas around M87, we restricted the 
analysis of the surface brightness to four different energy bands 
selected as follows: 

1 . Low energy band (0.4-0.9 keV) 

2. Fe-L line complex energy band (0.9-1.1 keV) 

3. Intermediate energy band (1.1-2.0 keV) 

4. High energy band (3.0-7.5 keV) 
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The principle is to compare the trend of fluxes in each pixel in 
each of the four bands to fluxes expected in these bands from 
model-spectra for a range of temperatures T={0.8, 1.0, 1.5, 1.7, 
2.0, 2.2, 2.5, 2.7, 3.0, 3.3, 3.6, 4.0} keV and metallicities z={0.1, 
0.3, .5, 0.7, 1.0, 1.3) solar (usi ng the relative abundance val- 
ues o flAnders & Grevessd(fl9 89)). Subsequently, the best-fitting 
model was determined using a least-chi-square fit. In all mod- 
els, the absorbing column density njj was set to 1.5 xlO 20 cm" 2 
dLieuet alJll996l) . and the redshift to 0.00436. 

We merged the two MOS detectors, which have very simi- 
lar spectral properties, and produced background-corrected flux 
maps in each of the four selected bands for the PN and com- 
bined MOS detectors. Each map was adaptively binned using a 
binning scheme based on the combined raw counts image from 
the three detectors in the 0.5-7.5 keV band. We performed one 
fit with a signal to noise ratio (SNR) of 50 and one with a SNR 
of 100 which ensured a minimum number of counts of 25... 100 
per bin in each of the four considered bands and especially the 
high energy 3.0-7.5 keV band which has the least counts. 

This method allowed us to clearly identify the eastern 
and southwestern arms of M87 as lower temperature features. 
Furthermore, a lower temperature component could be detected 
in the direction of the M87 jet (north-west of the core), almost 
in continuation of each of the arms, spatially connecting them. 
However, beyond a radius of around 11 arcmin, both the SN 
50 and SN 100 temperature maps became very noisy. This was 
likely due to the fact that our chosen energy bands were tuned to 
determine rather low temperatures accurately; for a good deter- 
mination of higher temperatures it would be necessary to sample 
the higher energy regime (above 3.0 keV) with more bands. The 
abundance pattern was also rather noisy throughout the produced 
maps; it seemed to be influenced significantly by the PN gap pat- 
tern, and did not show a clear distribution relating to the X-ray 
arms. This was likely due to the fact that our energy bands are 
very wide, leading to a washing-out of any emission lines with 
the exception of the Fe-L complex, which is not sufficient for an 
accurate abundance determination. 

In summary, this method was a very fast, computationally 
cheap way to gain insight into the spatial distribution of temper- 
atures in the observation. However, the result seemed to be rather 
sensitive to the choice of the energy bands, and the insufficient 
sampling of the higher-energy regime caused a noisy fit in the 
hotter regions of the gas halo. Moreover, it was an insensitive 
method for obtaining metallicity maps. 

2.2. Spectral analysis 

A more reliable method for creating temperature maps is a com- 
prehensive fit to the spectra from each bin using the XSPEC 
code, taking into consideration a large number of energy chan- 
nels and allowing the temperature and metallicity to vary con- 
tinuously until a best fit is found. But, given the large number 
of spatial bins for this observation, this method can be computa- 
tionally very demanding. 

Point sources were found using the source detection algo- 
rithm implemented in SAS and excised from the observation 
after a visual check to eliminate spurious detections. The cen- 
tral AGN was not considered a point-source contamination and 
was thus kept in the dataset throughout the analysis. We then 
adaptively binned the combined counts image in the 0.5-7.5 keV 
range using the weighed Voronoi tessellations method to a target 
signal-to-noise ratio of 100 and extracted the spectra of each bin 
defined in this way. The reason for binning the raw counts map 
and not the flux map is to obtain larger bins around the gap ar- 



eas where fewer photons are collected and to ensure an approxi- 
mately homogeneous accuracy in the temperature determination 
across the image. A SNR of 100 should give a temperature ac- 
curacy better than 5% (for a one-temperature fit) and a metal 
abundance accuracy of 10-15%. 

For the spectral ana lysis we chose the method described e.g. 
bv lArnaud et al.l (1200 ll) in which the vignetting correction is ac- 
counted for by adding a weight column to the event file, allowing 
the use of a single on-axis ancillary response file (ARF) contain- 
ing information about the telesope-specific effective area, quan- 
tum efficiency attenuation and filter transmission. The redistri- 
bution matrix files on the other hand were computed separately 
for each bin and each detector. Detector gaps were accounted for 
by dividing the flux in each bin through the corresponding ratio 
between the average unvignetted exposure time in the bin and 
the maximum unvignetted exposure time of the detector. 

The spectra in each bin were fitted in XSPEC using a 
MEKAL one-temperature model. The spectra from all the three 
detectors were grouped together and fitted simultaneously with 
the same best-fit model in order to have good statistics. We used 
the 0.5-7.5 keV range to determine the best fit. However, in the 
rare bins (approximately 10 bins out of a total of about 1500) 
where one of the EPIC detectors had an effective average ex- 
posure time of below 10% of the maximum, the data from that 
detector was omitted since the number of counts was too low to 
provide any reliable spectral information. This correction is nec- 
essary, as the gaps are not included in the vignetting correction 
information stored in the weight column. 

This method confirms the results obtained from the broad 
band fitting in that we can clearly identify the same substruc- 
ture as noted before: the eastern and southwestern arms of M87 
as lower temperature features, the lower temperature compo- 
nent NW of the nucleus in the direction of the M87 jet, and 
the nucleus of M87 distinguished as a high-temperature com- 
ponent. We also observe the positive temperature gradient in the 
regio n outside the arms, i n agreement with the radial analysis 
of Matsushita et al. (2002). The outer parts of the gas halo are 
significantly less noisy than derived from the broad band fit- 
ting. However, at smaller radii the values of the temperatures 
determined with the two methods agree very well. Our temper- 
ature map is also in good agreement with the temperature map 
obtained from the previous XMM- Newton observat i on of M 87, 
computed f ollowing the method of IChurazov et al.l (119961) and 
presented in lForman et al.l (12005b . 

The abundance pattern is much better defined with this 
method compared to the broad-band fitting, however some 
noise can still be seen. We note that the statistics in each 
bin were not sufficient to determine individual element abun- 
dances, therefore the abund ance ratios were assumed to be solar 
(lAnders & Grev esse 198 9)). Detailed ra dial profiles of the indi- 
vidual elements jMatsushita et al.ll2003l) however show that this 
is not generally the case. Given our signal-to-noise ratios, the 
obtained abundance values are most likely driven mainly by the 
iron content whose L-line complex is the most prominent fea- 
ture in the spectrum. We find that the abundance in the arms 
is in general lower than in the surrounding medium (Figure [3j, 
which is most likely due to the fact that we are using a single- 
temperature fit in a region which has been sho wn to have at 
least two components wi t h diffe r ent temperature s dBelsole et al.l 
(120011) . iMatsushita etail (120021) . iMolendil d2002l) ). As these au- 
thors have already noted, a one-temperature fit severely under- 
estimates the metallicity in this case. Moreover, by creating a 
map (not shown) of the absorbing hydrogen column density (ne) 
which was left free in the fit, we find elevated values in the 
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Fig. 3. Abundance map obtained from spectral analysis, using a 
Fig. 2. Temperature map obtained from spectral analysis, using binnins to SNR of 100 
a binning to SNR of 100. 



regions corresponding to the arms, which may also show that 
a single-temperature fit is not sufficient to appropriately model 
the data. We performed a second fit fixing the absorbing hydro- 
gen column density to 2.0xl0 20 cirT 2 which did not significantly 
change the temperature map. Therefore, any uncertainty in the 
ne seems to have little effect on the temperature results; in the 
following we will base our discussion on the initial temperature 
results obtained with an unconstrained ne in the fit. 



2.3. Pressure and entropy maps 

We used the temperature and spectrum normalization in each bin 
to determine a quasi-deprojected measure of the pressure and en- 
tropy, which we calculated as p — n e kT and S = kTn e m . The 
entropy as defined here is somewhat different from the thermo- 
dynamic quantity defined as S ~ (3/2)kln(Tp~ 2 ^) but is nev- 
ertheless a measure of adiabaticity and has become commonly 
used in cluster astrophysics. 

The electron density n e is determined from the spectrum 
normalization, which is defined and implemented in XSPEC as 
K = lO- H /[4nD 2 A (l+z) 2 ] J n 2 e dV, from which follows 



//, cv- ^l—D A (l +z) 

where D A is the angular diameter distance and z is the redshift. 
As can be seen from the above formula, we need to also calculate 
the volume V along the line of sight corresponding to each bin. 
This was determined as V ~ (4/3)D^Q(6» 2 u/ - # 2 ) I/2 , where £2 is 
the solid angle subtended by the bin and 8j„, 6 out are the angles 
corresponding to the smallest and respectively larges t distances 
between any of the bin pixels and the M87 nucleus dHenrv et al.l 



12004 . Thus, this method takes into account an approximate esti- 
mation of the three-dimensional extent of each bin, but assumes 
a constant temperature along the line of sight, and can therefore 
be viewed as a quasi deprojection. Since most of the emission 
in each bin originates from the densest gas which is found at 
the smallest effective 3D radii, the maps can be interpreted as 
measures of the pressure and entropy in a two dimensional slice 
through the middle of the cluster, perpendicular to the line of 
sight. 

We fitted the radial pressure and entropy profiles using non- 
parametric, locally weighted, linear regression smoothing. The 
data points and the corresponding smooth radial model are over- 
plotted in Figure [4] The pressure and entropy maps were di- 
vided by the resulting radial model in order to reveal small-scale 
fluctuations (Figures [5] and [6t. We ca n identify several features 
pointed out by iForman et alJ d2005h . such as the 3' (14 kpc) 
ring of enhanced emission corresponding to an increase in the 
pressure. This confirms the association of this substructure with 
a shock, while the NW/SE br ightness enhancements seen with 
Chandra (Formanetal. 2005) appear as low-entropy features. 
The X-ray arms are clearly low-entropy structures, which also 
seems to confirm the hypothesis that the gas in these regions was 
uplifted from more central regions where we expect the smallest 
entropy in hydrostatic equilibrium conditions. The next section 
presents a more detailed analysis of these features. 



3. Substructure in the M87 gas halo 

3.1. Temperature and Entropy 

The temperature and entropy maps both show similar substruc- 
ture details, so that we chose to discuss them in parallel. By far 
the most striking features in both maps are the E and SW arms 
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Fig. 5. Entropy deviations from a smooth radially symmetric 
model. The E and SW arms, as well as the low-entropy feature 
close to the NW of the core are seen in dark blue. The entropy 
edge to the S W is marked by the sharp green to yellow transition. 
A NW/SE asymmetry is easily seen in the map. 



Fig. 6. Pressure deviations from a smooth radially symmetric 
model. A relative pressure increase towards the SE and NW sug- 
gests an overall ellipticity in the pressure. A pressure decrease is 
found in the direction of the E and SW arms, which rise perpen- 
dicular to the SE/NW ellipticity long axis. A ring of enhanced 
pressure with a radius of roughly 3' is also seen. 



characterized by lower temperature and lower entropy with re- 
spect to the surroundings. The description of the thermal struc- 
ture of the arms will be the subject of an upcoming paper in 
which a more detailed spectral analysis will be employed in or- 
der to understand the complex multi-temperature stru cture that 
severa l authors have already f o und in these region s dMolendil 
d2002l) . lMatsushita et all d2002l) . lBelsole et all ( 12001b '). We point 
out that in the temperature and entropy maps both arms are seen 
to curve clockwise, thereby connecting to the larger-scale radio 
halos north and south of the nucl eus. This wa s alrea dy known 
for the SW arm (see for example iForman et alj ([2005)), but had 
not previously been observed for the eastern arm. 

Apart from the E and SW arms, another low-temperature 
and low-entropy feature can be identified close to the nucleus 
towards the NW, almost at the intersection of the lines defin- 
ing the general directions of the arms. This region corresponds 
to high X-ray luminosity and low entropy, implying a locally 
high gas densit y. Additionally, enh anced Ha emission has been 
observed here (Spa rks et al.ll2004b . which makes this region an 
attractive candidate for hosting a classical cooling flow and as- 
sociated mass deposition. However, detailed spectral analysis of 
the region, as described in the next section, does not provide evi- 
dence for cooling to the lowest X-ray detectable temperatures at 
the rates predicted in the absence of any heating sources. 

Another feature found in both maps is an edge to the SE at 
a radius of roughly 6'. As one moves outwards across this edge, 
both the temperature and entropy increase, therefore the proper- 
ties at the jump are more consisten t with the cold-front interpre- 
tation of iMarkevitch et a l. (2001) than with the characteristics 
of a shock. A more detailed discussion of this is presented in the 
next section. 



Finally, a feature easily seen in the entropy map but which is 
not readily apparent in the temperature map is a NW/SE asym- 
metry. Within a radius of 6' the entropy values to the NW are 
clearly more elevated than in the SE. The edge to the NW where 
the entropy decreases corresponds spatially very well to the edge 
of the NW large-scale radio bubble. Possible physical explana- 
tions for this are also presented in the next section. 

3.2. Pressure 

A map of the pressure deviations from a radially symmetric, 
smooth model shows three noteworthy features. The first of these 
is a relative pressure increase towards the NW and SE. Since we 
subtracted a radially symmetric model, this would suggest an 
overall ellipticity of the pressure distribution and consequently 
of the underlying dark matter distribution under the assumption 
of hydrostatic equilibrium. Secondly, in the direction of the E 
and SW X-ray arms (coinciding with the prominent radio lobes 
and orthogonal to the regions of enhanced pressure), we find a 
pressure decrease. Thirdly, there is a ring of enhanced pressure 
with a radius of roughly 3', coinciding with the position of the 
weak shock proposed bv lForman et alj d2005h . 

We tried to confirm and quantify the ellipticity in the initial 
pressure map by fitting to it elliptical isocontours with the "el- 
lipse" task in IRAF. The position angle and ellipticity were in- 
dependently fitted for each pressure contour and the model need 
not be smooth in intensity as a function of the semi-major axis. 
We masked the regions corresponding to the radio lobes from the 
fit, since here it is likely that the pressure deviates from an ellip- 
tical model. Dividing the initial pressure map by the fitted model 
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Fig. 7. Pressure deviations from an elliptical model (ratio). 
Contour levels are drawn at 0.9, 0.95, 1.05 and 1.1. Two regions 
roughly corresponding to the end of the E arm and to the SW 
arm are seen to have lower than average values, indicating the 
possible presence of nonthermal pressure support. 



obtained with IRAF, we are able to identify deviations from an 
elliptical rather than a radially symmetric model (Figure |7J. We 
find that, with respect to the elliptical model, all substructure 
seen previously no longer appear, with the exception of the lower 
pressure in regions that coincide with the radio lobes. Here the 
decrease in observed pressure with respect to the elliptical model 
is on average 5-10%. Assuming an overall pressure balance, this 
decrease suggests that in the radio lobes there should be an addi- 
tional contribution to the pressure. This pressure source may be 
provided by relativistic electrons injected by the AGN. The ring 
at 3' is no longer visible since we did not require the pressure 
intensity in the model to be smooth with increasing semi-major 
axis, while the higher pressures NW and SE of the nucleus which 
appeared as deviations from a radially symmetric model can be 
entirely explained by an elliptical pressure distribution. 

4. Discussion 

4.1. Absence of a classical cooling-flow on smaller scale 

While a large-scale cooling flow has been ruled out, cooling and 
condensation could still happen locally, since it is not easy to 
exactly balance heating and cooling everywhere throughout the 
cooling core region. A prime candidate target to look for such a 
local cooling flow region is the small low-entropy feature NW 
of the nucleus where also diffuse optical emission lines are ob- 
served. 

We extracted a spectrum from this region and fitted it in 
XSPEC with a vmcflow model fixing the low temperature cutoff 
at the minimum level available in the model (8 1 eV) in order to 
probe the existance of a classical cooling flow spectrum. We find 
a very poor fit, especially around the Fe-L complex where the 
data lie below the spectral model between 0.6-1 keV (Figure[8]l. 
A vmekal+vmcflow model with a fixed low-temperature cutoff 
at 8 1 e V provides a good fit for an ambient temperature (vmekal) 
of 1.508 + 0.027 keV and a mass deposition rate of 3.18 + 0.16 




2 5 
channel energy (keV) 



Fig. 8. Data and best-fit vmcflow model with fixed low- 
temperature cutoff for the small-scale low-temperature feature 
NW of the core. The disagreement of the model with the data is 
easily visible, indicating the absence of a classical cooling flow. 




1 1 2 5 

channel energy (keV) 

Fig. 9. Data and best-fit two-temperature vmekal model for the 
small-scale low-temperature feature NW of the core. This model 
is clearly in better agreement with the data than the classical 
cooling flow model. 

xl0~ 2 Mo/year. Accordi ng to a classical cooling flow analysis 
done by Bohringer ( 1999), the mass deposition rate in this region 
in the absence of a heating source should be 0.1-0.2 M /year, 
which is a factor of about 5 higher than what we find from spec- 
tral fitting. A two-temperature vmekal model with Tl = 1.016 
± 0.022 and T2 = 1.855 + 0.030 also describes the data very 
well. Therefore, although the density in this smaller-scale region 
is very high and although the region is associated with Ha emis- 
sion, we find that it does not exhibit spectral evidence of a clas- 
sical cooling flow and sufficient mass deposition rate. 

4.2. Motion of the galaxy to the NW 

To investigate the NW/SE asymmetry seen in the entropy map, 
we extracted spectra from concentric ring sectors towards the 
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NW (-25 to 135 degrees counterclockwise of West) and SE (200 
to 280 degrees counterclockwise of West). The angles are chosen 
such that the E and SW arms are excluded. The spectra were fit- 
ted with a on e-temperature vmekal mo del, since we expect from 
the results of iMatsushita et al.l d2002l) and iMolendil d2002|) that 
the gas outside the arm regions is single-phase locally. The re- 
sults plotted in Figure[TT]show that the radial temperature profile 
towards the SE lies below the profile towards the NW, while the 
radial abundance profile reveals higher values in the SE than in 
the NW. The simplest phenomenon that explains such a NW/SE 
temperature and abundance asymmetry is a motion of M87 to the 
NW through the hot gas halo, or alternatively, a large-scale bulk 
motion of the hot gas halo towards the SE while the galaxy itself 
remains in place. In either case, there are two possible mecha- 
nisms related to the motion which would generate the observed 
asymmetry. 

The first possibility is that colder, highly-abundant gas near 
the center is ram-pressure stripped during the motion. In the case 
of a large-scale bulk motion of the gas halo, the stripped gas 
is then carried towards the SE along with the bulk flow, while 
if the galaxy is moving towards the NW it is left behind in a 
wake. However, the pressure map is symmetrical and shows no 
increase to the NW as one would expect if ram-pressure were 
important. 

The second possibility, which is favoured by our entropy and 
pressure maps, is that, due to the relative velocity between the 
nucleus and the gas halo, bubbles emitted initially to the north- 
west are advected downstream. This effect is important if the 
relative bulk velocity is at least as high as, and of the order of, 
the bubble rising velocity and would imply that a large part of the 
bubbles emitted by the AGN eventually rise to the SE indepen- 
dently of their initial direction of emission. Much more mixing 
is therefore induced towards the SE, favoring the transport of 
metals, which explains the relative iron and silicon abundance 
increas e in this direction (Fig ure [TTI (b) and (c)). Recent simu- 
lations Roedig er et al.l (H006) indeed show that bubble-induced 
metal transport results in elongated abundance profiles along the 
direction of propagation of the bubbles. The observed NW/SE 
temperature asymmetry can be explained using this bubble- 
advection scenario if we consider that the bubble-induced mix- 
ing not only favors the transport of metals in the atmosphere but 
can also bring colder gas from the center out to larger radii. It has 
been already noted (e.g. Nul sen et al.l d2002l) ) that radio bubbles 
are often surrounded by rims of cold, dense gas which they uplift 
during their inflation and rising. Therefore it is probable that a 
preferential propagation of a series of bubbles to one direction 
may result in lowering the downstream average temperature. 

Further indication that bu bbles are advected d ownstream is 
provided by Chandra data in iForman et al.l d2006l) : the authors 
describe a series of four consecutive bubbles to the SSE while 
no counterpart of such structure is seen in the opposite direction. 

Finally, one other observational feature leads us to believe 
that M87 is moving to the NW. If we consider the analogy of 
M87 with a typical wide-angle tailed (WAT) radio source, the 
motion through the intracluster medium may help to explain why 
the E and SW radio lobes, and the corresponding X-ray arms, do 
not lie along directions exactly opposite to each other. However, 
the wide angle between the lobes and especially the radio "ear" 
seen at the end of the eastern arm suggest that the buoyant veloc- 
ity is at least comparable to the bulk flow velocity. This agrees 
qualitatively with our previous statement that the relative bulk 
velocity should be on the order of the bubble rising velocity. 

Let us recons ider the mode l of ICM heating through radio-jet 
inflated bubbles. iBirzan et alj d2004l) . among others, show that 
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Fig. 11. Temperature and abundance NW/SE asymmetry. The 
solid line represents the profile NW of the core, the dashed red 
line SE of the core. 



the bubble enthalpy is in most cases e nough to halt the c ool- 
ing flow. However, the picture drawn bv lNulsen et al.l d2002l) for 
Hydra A and in this work also for M87 is much more com- 
plex. While the AGN does inject energy into the ICM in the 
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form of bubble enthalpy, bubble formation is likely often as- 
sociated with the uplift and mixing of cool gas from the cen- 
ter. A well known example for the uplift are the X-ray arms of 
M87 (IChurazov et al.ll2001h . In this paper we also find for the 
first time evidence of this phenomenon at smaller temperature 
and density contrasts with respect to the ambient medium within 
the same galaxy, contrasts which could therefore not have been 
observed without the high-quality statistics of the dataset. This 
suggests that cool-gas uplifting by buoyant bubbles may be more 
wide-spread than previously thought, and has not been detected 
so far solely owing to insufficient statistics. Whether the uplifted 
cool negatively buoyant gas will fall b ack towards the cl uster 
center thermalizing its potential energy dNulsen et al.ll2002h . or 
whether it is heated by thermal conduction once it reaches larger 
radii, is a matter for future models; however, it seems that bub- 
bles do generally transport cool gas out of the center. 

4.3. Possible core oscillations 

While the abundance and temperature profiles of the X-ray gas 
relatively close to M87 (within the inner ~ 6') are well explained 
by a simple relative motion between the gas and the galaxy, a 
look at the features seen further towards the outskirts reveals a 
more complex nature of the ongoing physical phenomenon. The 
model-divided entropy map reveals an edge at about 6' SE of the 
core which is suggestive of a cold front, with the pressure across 
the front staying roughly constant and the temperature and en- 
tropy being lower on the more X-ray luminous side. This feature 
corresponds in the metallicity map (Figure [3) to a visible sharp 
edge in the abundance distribution, emphasizing the presence of 
a contact discontinuity here. Another relatively sharp edge in the 
metallicity can also be seen at about 3' NW of the core, although 
it is not associated with any visible entropy or temperature fea- 
tures. Moreover, also NW of the core, we find a region of lower 
entropy beyond 6'. 

The opposite and staggered placement of these low-entropy 
regions and abundance edges suggests that the current relative 
motion of the gas and the galaxy, proposed in the previous sec- 
tion, may be part of a lon ger process of oscillat ory motions 
along the NW/SE direction. Churazov et al.l (120031) observed a 
similar E/W asymmetry and a set of sharp edges in entropy, 
temperature and surface brightness in the Perseus cluster and 
performed simulations which suggest that minor mergers may 
trigger oscillations of the cluster gas, thus generating multi- 
ple sharp edges on opposite sides of the central galaxy along 
the direction of the merger. Similar opposite and staggered fea- 
tures generated by oscillations of the intracluster gas are seen 
also in numerical simulat i ons bylTittlev & Henriksenl (120051) and 
lAscasibar & Markevitchl d2006h . 

Further evidence of a possible oscillatory motion of the gas 
in the M87 halo comes from an analysis of the pressure map. 
Fitting an elliptical model to the raw pressure map as described 
in Section 4.5 and allowing the centers of the elliptical isobars 
to vary, we find that these centers systematically shift towards 
the SE for higher values of the semimajor axis. Thus the galaxy 
potential located at smaller scales is not exactly centered with 
the large-scale halo potential, which would lead to the oscilla- 
tion suggested here. In a more detailed future analysis, this can 
be used to reveal information about the amplitude, energy, and 
time-scale of these oscillations, which may help in understand- 
ing how this phenomenon contributes over time to the heating 
vs. cooling balance, and how it competes with other proposed 
methods for energy input into the ICM. 
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Fig. 12. Entropy and temperature radial profiles, calculated as an 
average of the values in the bins whose center falls inside circular 
annuli sectors with an opening angle of 90 degrees NW and SE 
of the core. The errorbars represent the corresponding standard 
deviation in each annulus sector. 



4.4. Large-scale correlation of the radio emission and 
entropy map 

One of the striking large-scale features in Figure [10] is the spa- 
tial coincidence of the northern outermost radio contour-edge 
and a decrease in the entropy by 10-15%, which cannot be ex- 
plained either by a directe d motion of M87 t o the NW, or by the 
possible core oscillations. IOwen~et al. (2000) note that the radio 
contour-edge represents a sharp decrease in radio intensity over 
a range of wavelengths, therefore indicating interaction between 
the radio-loud plasma and the hot gas halo. 

To illustrate the entropy variation across the radio edge, we 
averaged the entropy in circular annuli sectors with an opening 
angle of 90 degrees in the NW and SE quadrant respectively. We 
compare the radial profiles of the entropy in these two quadrants 
in Figure [12] which also shows radial trends of the temperature 
in the corresponding regions. Up to a radius of roughly 6'-7' the 
entropy in the NW quadrant is higher than in the SE quadrant. 
Beyond this radius, which corresponds roughly to the radio edge, 
the entropy and temperature NW of the core show a much shal- 
lower increase with radius than in the SE, the temperature dis- 
playing even a small negative gradient. The spatial coincidence 
between the north radio edge and the point where the entropy and 
temperature gradients NW of the core decrease suggests a link 
between the injection of radio-loud plasma into the ICM and the 
processes leading to entropy generation and gas heating. Such a 
link is at the basis of the feedback models in which mechanical 
input, usually by a radio source, offsets the cooling flow, heat- 
ing the gas at the cluster center and flattening the central entropy 
profiles. The exact details of the entropy-generation mechanism 
remain unclear. 
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4.5. Ellipticity of the underlying dark matter potential 



As mentioned above, we fitted the raw pressure map with an 
elliptical model in IRAF in order to quantify the eccentricity 
and position angle-dependency with semi-major axis. The re- 
sults of the obtained fit show larger variations in the central part 
but outside 3' a rather constant position angle and ellipticity are 
reached. The position angle values outside 3' range between 140 
and 150 degrees, measured counterclockwise of north. These 
values are somewhat smaller than, but compa rable to, results 
from t he analysis of an optical image of M87 bv lCarter & Dixonl 
( 1978), who obtained values for the elliptical isophote-fit of 160 
+ 5 deg rees for the position an gle. Our results agree also with the 
work o flBohringer et all (1 1997b on ROSAT PSPC data, giving an 
ellipticity of the X-ray surface brightness with a position angle 
of 158 + 10 degrees. Outside of 3', the ellipticity varies generally 
between 0.08 and 0.17, higher values of 0.13-0.17 being found 
in the regions from which the radio lobes h ad been cut out. This 
is in agreement with Bohring er et al.l (Q997), who found an ellip- 
ticity of 0.1 to 0.16. The ellipticity in the X-rays thus seems to be 
significantly smaller than in the optical regime, where it ranges 
from 0.3-0.5. This can be explained by the fact that anisotropic 
velocity dispersions of the stars can sustain a higher ellipticity 
of the galaxy in the optical domain while the isotropic nature of 
the X-ray emitting gas halo makes it more round. Surprisingly, 
despite the evidence of core oscillations seen in the entropy map, 
the position angles of the X-ray isobars and the optical isophotes 
agree very well, and the optical ellipticity is as expected roughly 
3 time s higher than the X-ray value (see also lBuote & Cani zares 
(l!996h V Therefore the deviations of the hot gas atmosphere from 
hydrostatic equilibrium cannot be very large outside the arm re- 
gions. 

A further look at the position angle values found in the fit 
reveals the fact that the radio lobes rise in directions orthogo- 
nal to the semimajor axis of the elliptical dark matter profile, 
thereby following the steepest dark-matter potential gradient as 
one would expect. 

4.6. Heating of the ICM by weak shocks 

As seen in Figure [6] we are able to sp e ctroscopically confir m the 
3' shock seen bv lForman eTail (l2005l) . lForman et al.l (f2006) with 
Chandra. Weak shocks are becoming widely accepted as candi- 
dates for heating mechanisms of the X-ray gas, therefore a de- 
tailed understanding of their properties is required. To describe 
the M87 3' shock more quantitatively, we used the elliptical iso- 
bars from the fit described above to define eight spectral extrac- 
tion regions with semi-major axes between 2 and 4 arcminutes. 
The temperature and spectrum normalization corresponding to 
the regions contained between each two consecutive elliptical 
isobars in the chosen semi-major axis range were determined 
using a vmekal one-temperature spectral fit. The corresponding 
pressure was then determined for each region as described in 
Section 2.3. The results are plotted in FigureQj] In the plot, the 
pressure was divided by a smooth model to emphasize the jump. 

The temperature shows a jump of ~0.1 keV, which relative 
to the overall temperature of ~2 keV, represents ~5%. Using 
the Rankine-Hugoniot jump conditions assuming y = 5/3, this 
amounts to a Mach number of 1.05. The pressure jump amounts 
to 8-10% with respect to the smooth model, which implies a 
Mach number of about 1 .04, in agreement with the Mach num- 
ber we derived from the temperature jump. This confirms the 
supersonic nature of the weak shock at 3', and seems consis- 
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Fig. 13. Projected temperature and relative pressure jumps 
around the 3' shock, corresponding to Mach numbers of 1.05 
and 1 .04 respectively. 



tent with the Mach numbe r M«; 1 .2 calculated f rom deprojection 
analysis of Chandra data (Formanet al. 2006) considering that 
the Mach numbers derived from this analysis are lower limits to 
the real values due to projection effects and the smoothing effect 
of the XMM PSF Accurate spectral deprojection analysis is the 
scope of upcoming work. 



5. Summary and conclusions 

We present the results from a 109 ksec XMM-Newton observa- 
tion of the hot gas halo surrounding the giant elliptical galaxy 
M87. We use two methods, namely broad-band fitting and full 
spectroscopy, to create spatially resolved temperature maps from 
the data, and present the advantages and disadvantages of each 
method. The results of the spectroscopic analysis are used to pro- 
duce entropy and pressure maps. We describe and discuss the 
features seen in these maps, the most important of which being 
the cool E and SW X-ray arms, which coincide with powerful 
radio lobes, a weak shock ring with a radius of 3', an overall el- 
lipticity in the pressure map and a NW/SE asymmetry in the en- 
tropy map. For the 3' weak shock we find jumps in both the pres- 
sure and temperature corresponding to Mach numbers of 1.04 
and 1.05 respectively, which we interpret as lower limits to the 
real Mach number values. The pressure map shows an overall el- 
lipticity of between 0.08 and 0.17, with position angles in good 
agreement with the optical data. Under the assumption of hy- 
drostatic equilibrium, the ellipticity in the pressure map implies 
an elliptical underlying dark matter potential. Furthermore, we 
find that the X-ray arms and associated radio lobes rise roughly 
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perpendicular to the semimajor axis of the pressure ellipticity 
following the steepest gradient of the dark matter potential. 

The NW/SE entropy asymmetry is indicative of the motion 
of M87 through the surrounding hot gas, associated with down- 
stream advection of the bubbles injected by the AGN. We con- 
clude that bubble-induced cold gas mixing from the center may 
appear on a wide variety of scales and should be considered into 
the cooling-flow energy balance calculations. NW and SE edges 
in the entropy map moreover suggest that the possible motion of 
M87 to the NW is part of a longer cycle of core oscillations. 

Acknowledgements. We would like to thank W. Forman, P. Nulsen, E. Churazov 
and G. W. Pratt for helpful discussion. We acknowledge the support by the 
DFG grant BR 2026/3 within the Priority Programme "Witnesses of Cosmic 
History". The XMM-Newton project is an ESA Science Mission with instru- 
ments and contributions directly funded by ESA Member States and the USA 
(NASA). The XMM-Newton project is supported by the Bundesministerium 
fuer Wirtschaft und Technologie/Deutsches Zentrum fuer Luft- und Raumfahrt 
(BMWI/DLR, FKZ 50 OX0001), the Max-Planck Society and the Heidenhain- 
Stiftung, and also by PPARC, CEA, CNES, and ASI. AF acknowledges support 
from BMBF/DLR under grant 50 OR 0207 and MPG. This project is partially 
supported by a NASA grant NNG05GH40G to UMBC. AS also thanks Steven 
Diehl and Thomas Statler for making their Weighted Voronoi binning algorithm 
publicly available. 



References 

Anders, E. & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 
53, 197 

Arnaud, M., Neumann, D. M., Aghanim, N., et al. 2001, A&A, 
365, L80 

Ascasibar, Y. & Markevitch, M. 2006, ApJ, 650, 102 
Belsole, E., Sauvageot, J. L., Bohringer, H., et al. 2001, A&A, 
365, LI 88 

Btrzan, L., Rafferty, D. A., McNamara, B. R., Wise, M. W, & 

Nulsen, P. E. J. 2004, ApJ, 607, 800 
Bohringer, H. 1999, in Diffuse Thermal and Relativistic 

Plasma in Galaxy Clusters, ed. H. Boehringer, L. Feretti, & 

P. Schuecker, 115-120 
Bohringer, H., Belsole, E., Kennea, J., et al. 2001, A&A, 365, 

L181 

Bohringer, H., Matsushita, K., Churazov, E., Ikebe, Y., & Chen, 

Y. 2002, A&A, 382, 804 
Bohringer, H., Neumann, D. M., Schindler, S., & Huchra, J. P. 

1997, ApJ, 485, 439 
Bruggen, M., Ruszkowski, M., & Hallman, E. 2005, ApJ, 630, 

740 

Buote, D. A. & Canizares, C. R. 1996, ApJ, 457, 565 
Burns, J. O. 1990, AJ, 99, 14 

Cappellari, M. & Copin, Y. 2003, MNRAS, 342, 345 

Carter, D. & Dixon, K. L. 1978, AJ, 83, 574 

Churazov, E., Bruggen, M., Kaiser, C. R., Bohringer, H., & 

Forman, W. 200 1 , ApJ, 554, 26 1 
Churazov, E., Forman, W., Jones, C, & Bohringer, H. 2003, 

ApJ, 590, 225 

Churazov, E., Gilfanov, M., Forman, W., & Jones, C. 1996, ApJ, 
471,673 

Cowie, L. L. & Binney, J. 1977, ApJ, 215, 723 
David, L. P., Nulsen, P. E. J., McNamara, B. R., et al. 2001, ApJ, 
557, 546 

Diehl, S. & Statler, T. S. 2006, MNRAS, 368, 497 
Fabian, A. C. & Nulsen, P. E. J. 1977, MNRAS, 180, 479 
Fabian, A. C, Sanders, J. S., Allen, S. W., et al. 2003, MNRAS, 
344, L43 

Forman, W., Churazov, E., Jones, C, et al. 2006, ArXiv 
Astrophysics e-prints [astro-ph/0604583] 



Forman, W., Nulsen, P., Heinz, S., et al. 2005, ApJ, 635, 894 
Harms, R. J., Ford, H. C, Tsvetanov, Z. I., et al. 1994, ApJ, 435, 
L35 

Heinz, S., Brueggen, M., Young, A., & Levesque, E. 2006, 

ArXiv Astrophysics e-prints [astro-ph/0606664] 
Henry, J. P., Finoguenov, A., & Briel, U. G. 2004, ApJ, 615, 181 
Lea, S. M., Silk, J., Kellogg, E., & Murray, S. 1973, ApJ, 184, 
L105 

Lieu, R., Mittaz, J. P. D., Bowyer, S., et al. 1996, ApJ, 458, L5 
Markevitch, M., Vikhlinin, A., & Mazzotta, P. 2001, ApJ, 562, 
L153 

Matsushita, K, Belsole, E., Finoguenov, A., & Bohringer, H. 

2002, A&A, 386, 77 
Matsushita, K, Finoguenov, A., & Bohringer, H. 2003, A&A, 

401,443 

McNamara, B. R. & O'Connell, R. W. 1989, AJ, 98, 2018 
Molendi, S. 2002, ApJ, 580, 815 

Nulsen, P. E. J., David, L. P., McNamara, B. R., et al. 2002, ApJ, 
568, 163 

Owen, F. N., Eilek, J. A., & Kassim, N. E. 2000, ApJ, 543, 611 
Peterson, J. R., Kahn, S. M., Paerels, F. B. S., et al. 2003, ApJ, 
590, 207 

Peterson, J. R., Paerels, F. B. S., Kaastra, J. S., et al. 2001, A&A, 
365, L104 

Pratt, G. W. & Arnaud, M. 2002, A&A, 394, 375 

Read, A. M. & Ponman, T. J. 2003, A&A, 409, 395 

Roediger, E., Bruggen, M., Bohringer, H., Rebusco, P., & 

Churazov, E. 2006, MNRAS submitted 
Sparks, W. B., Donahue, M., Jordan, A., Ferrarese, L., & Cote, 

P. 2004, ApJ, 607, 294 
Stewart, G. C, Fabian, A. C, Jones, C, & Forman, W. 1984, 

ApJ, 285, 1 

Tittley, E. R. & Henriksen, M. 2005, ApJ, 618, 227 
Tonry, J. L., Dressier, A., Blakeslee, J. P., et al. 2001, ApJ, 546, 
681 

Watson, M. G, Augueres, J.-L., Ballet, J., et al. 2001, A&A, 
365, L51 

Werner, N., Bohringer, H., Kaastra, J. S., et al. 2006, A&A, 459, 
353 



A. Simionescu et al.: The gaseous atmosphere of M87 seen with XMM-Newton 



11 



200 



L 
u 

> 

0) 



O 

c 

CD 




20 40 60 
radius (kpc) 



SO 



< 

u 



X) 



3 
CO 
CO 
■!■ 




00.0 



(a) Entropy profile and fitted smooth, non-parametric radial model (b) Pressure profile and fitted smooth, non-parametric radial model 

Fig. 4. Data points and corresponding smooth radial model by which the entropy and pressure maps were divided in order to reveal 
small deviations from spherical symmetry 




Fig. 10. Entropy (left) and pressure (right) deviations from a smooth radially symmetric model overlaid with radio contours (90 cm). 
Radio map kindly provided by F. Owen. The E and SW radiolobes coincide well with the regions of low entropy in the X-ray arms. 
Also, the edge of the large radiolobe to the north roughly coincides with a NW edge in the entropy map. 



